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ABSTRACT 

Recent  work  has  shown  how  to  test  a  simple  hypothesis 
of  uniformity  on  the  interval  (0,1)  by  using  estimates  of 
entropy.  In  this  paper  we  use  Monte  Carlo  methods  to 
extend  previous  tables  of  critical  points  and  power  for 
such  entropy  tests  to  the  large  sample  sizes  likely  to  ;e 
desirable  when  evaluating  the  output  of  one  or  more 
random  number  generators.  A  comparison  with  asymptotic 
critical  points  and  power  is  made.  The  results  are  used 
to  evaluate  a  number  of  commonly  used  random  number 
generators.  At  least  one  random  number  generator  judged 
"good"  by  the  spectral  test  (multiplier  5**15)  is  found 
unsuitable  for  use. 
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1. 


INTRODUCTION 


An  estimator  of  entropy  was  proposed  by  Vasicek 
(1976)  and  studied  by  Dudewicz  and  van  der  Meulen  (1979), 
who  showed  how  this  estimator  could  be  used  to  test 
whether  a  random  sample  comes  from  the  uniform  distribu¬ 
tion  on  (0,1)  and  gave  Monte  Carlo  and  asymptotic  evalua¬ 
tions  of  the  percentage  points  of  the  distribution  of  the 
estimator.  In  this  paper  we  utilize  this  test  for  uni¬ 
formity  to  evaluate  random  number  generators .  A  compari¬ 
son  of  asymptotic  with  Monte  Carlo  for  percent  points  and 
power  is  also  given. 

The  differential  entropy  (or  entropy)  of  a  random 

variable  X  with  density  function  f  is  defined  as 

00 

H(X)  =  -  /  f (x)  log  f(x)dx  (1.1) 

—  00 

and  has  the  properties  that  if  Xe[0,l]  w.p.  1  then 

H(X)  <  0,  (1.2) 

and  among  all  densities  f  concentrated  on  [0,1]  the 
uniform  fQ  maximizes  H(f)  to 

H( f q )  =  0.  (1.3) 

The  following  two  definitions  were  given  by  Dudewicz  and 
van  der  Meulen  (1979),  and  use  the  above. 
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Definition  (1.4).  Two  densities  and  f2  are  said  to  be 
entropy-distinguishable  if  H(f^)  4  Hffj). 

Definition  (1.5).  A  density  f*  is  said  to  be  entropy- 
unique  (or  e-unique )  in  a  class  _  of  densities  if 
f*  e  .  and  $  f  e  i  ,  f  4  f * ,  such  that  H(f)  =  H(f*). 


By  property  (1.3),  we  know  that  in  the  class  :  of 
densities  on  (0,1),  fg  is  e-unique.  This  fact  will  be 
central  below. 

Now  let  X ^ ,  ...  ,  Xn  be  a  random  sample  from  an 
absolutely  continuous  distribution  F  with  density  f  and 


let  X, 


be  the  order  statistics.  The 


‘(1)  ’  ‘  ’  X(n) 

estimator  of  H(f)  (Vasicek  (1976))  is 

n 


H 


m 


,n  =  n  1  J=1  log{75  CX(i+m)  '  X(i-m)]} 


(1.6) 


where  1  £  m  <  n/2,  =  X^-jj  for  j  <  1,  and  X(j)  =  ^(n) 

for  j  >  n . 

Dudewicz  and  van  der  Meulen  (1979)  proposed  a  test 
for  uniformity  based  on  the  e-uniqueness  of  fg.  We  are 
concerned  with  this  test  and  so  describe  it  here. 

Let  X^,  ...,  Xn  be  a  random  sample  from  an  absolutely 
continuous  distribution  F  with  density  f  concentrated  on 
[0,1].  Let  fg  denote  the  U(0,1)  density.  The  level  a  test 


rejects  H 


0 


f  =  f  in  favor  of  :  f  4  fg  if  and  only  if 


3 


H 


m 


,n  i  Ha*(m ,n > 


(1.7) 


where  Ha*(m,n)  is  the  100a  percentile  point  of  the  distri¬ 
bution  of  H  under  f. 


m  ,n  ---  -  -  o  ' 

By  the  e-uniqueness  of  among  all  densities  con¬ 
centrated  on  [0,1]  and  by  the  consistency  of  the  estimator 
(Vasicek  (1976)),  it  follows  that  the  above  test  is  con¬ 
sistent  against  all  alternatives  f  on  [0,1].  Further,  in 
Du.dewicz  and  van  der  Meulen  (1979)  it  is  shown  that  if 


f  is  concentrated  on  [0,1]  then,  w.p.  1,  H  <  0. 

m  ,n  — 


In  Dudewicz  and  van  der  Meulen  (1979),  the  above  test 
is  studied  in  detail  using  both  analytical  and  Monte  Carlo 
techniques.  The  latter  are  necessary  since  the  form  of 
the  distribution  of  Hm  n  appears  to  be  analytically 
intractable.  From  Dudewicz  and  van  der  Meulen  (1979)  we 
know 


Theorem  1.8.  Under  HQ  ,  if  m  =  o(n1/3-<5),  6  >  0,  the 

quantities 


<6mn)1/2[Hmin  -  log(^>  +  log(n  +  1)  +  Y  -  R(l,  2m-l)J  (1.9) 


'2m 


and 


,1/2, 


(6mn)  [H^n  +  log(2m)  +  y  -  R(l,  2m-l)] 


(1.10) 


are  asymptotically  N(0,1)  as  n 
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Here  y  is  Euler's  constant,  y  =  0.S772...  and,  for 
3  <  k,  R(j,  k)  =  k”1  +  (k-1)"1  +  ...  +  j"1.  Using  (1.9) 
we  then  have  that 

=  ^  ^a)(6mn)  -  log(2m)  -  y  +  R(l,  2m-l) 

-  log  (1.11) 

is  an  asymptotic  approximation  to  the  100a  percentile  point 
of  the  distribution  of  Hm  n  where  $(•)  is  the  standard 
normal  distribution  function. 

For  asymptotic  power,  Dudewicz  and  van  der  Meulen 
(1979)  show 


Theorem  1.12.  For  any  bounded  positive  step  function 
alternative  density  f^  on  (0,1), 

nl/2(hT7T2m)1/?[Hm,n  +  log(2m>  +  Y  -  R(l.  2m-l)  -  H(f1)] 

is  asymptotically  normal  with  mean  zero  and  variance 

x(2m)  +  Varf  log  f , ,  where 
X1  1 


it  /  6  -  1 


x  (k )  = 


if  k  =  1 


I 

1  (2k2-2k+l){I^-  -  (1  +  + 


(k-1) 


y}-2k+l 


if  k  >  2  . 


Hence  an  asymptotic  approximation  to  the  power  of  a 
level  a  test  under  alternative  f  is 


2. 


ONE-SAMPLE  EVALUATION  BASED  ON  ASYMPTOTIC  PERCENTILE 


RANKS 


We  now  use  test  (1.7)  to  evaluate  whether  each  of 
nine  of  the  uniform  random  number  generators  available 
in  Dudewicz  and  Ralley  (1981)  have  entropies  consistent 
with  uniformity.  At  the  same  time,  the  form  of  the 
test  (1.7)  suggests  that  if,  for  two  random  number 
generators  ...  ^  and  ~  ^  > 


H  (  ...  ) 
m,n  1 


H 


m  ,n 


then  is  "better  than"  .  This  provides  a  basis  for 

comparative  evaluation  of  the  random  number  generators. 

Thus  we  first  obtain  H  (.-/.)  for  each  of  nine  of  the 

m,n  i 

random  number  generators  h  in  Dudewicz  and  Ralley  (1981): 
KNCG  (using  IX  =  1 ,  LI  =  452807053,  C  =  0.0,  PIP  = 

2. **31),  SRAND ,  KERAND,  UNI,  RN1 ,  RN2 ,  RN3 ,  RN4  and  RN5 . 

Fixing  n  =  10,000  and  letting  m  =  1(1)10(5)20(10)40, 

comparison  between  the  results  obtained  by  Monte  Carlo 
methods  and  by  asymptotic  methods  is  one  of  our  goals. 

This  value  of  n  automatically  excluded  RN5  from  considera¬ 
tion  for  m  <  5  since  the  period  of  RN5  is  only  2048  (  =  2^), 
hence  leading  to  definite  rejection  due  to  gaps  of  size 
zero  for  small  m. 

For  each  of  the  0  generators,  the  statistics 
Hm  ^oooo  were  computed  from  a  single  sample  of  10,000 
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numbers  for  14  values  of  m  (10  for  RN5)  as  above.  The 
percentile  ranks  of  these  statistics,  expressed  in  terms 
of  their  asymptotic  distributions  under  (see  equation 
(1.10)),  were  plotted  in  Figure  1. 

Figure  1 

Examination  of  that  graph  reveals  that,  among  all  the 
generators  considered,  RN2  ,  RN4 ,  and  KERAND  have  entropy 
estimates  most  consistent,  almost  too  consistent,  with 
uniformity,  while  UNI  gives  good  entropy  estimates  in  the 
same  sense.  The  poorest  generator  is  RNCG,  while 
the  rest  performed  reasonably  well.  (To  be  rejected,  the 
n  -  value  of  a  generator  needs  to  have  a  percentile 
rank  5.00,  at  the  top  of  the  graph  scale.)  These 
interpretations  are  based  on  a  single  sample  of  size 
10,000  and  therefore  should  be  regarded  more  as  an 
illustration  of  use  of  n  for  comparing  random  number 
generators  than  as  a  definitive  statement  regarding  them. 
More  replications  are  used  below  in  a  way  reminiscent  of 
the  extensive  testing  of  random  number  generators  reported 
in  Dudewicz  and  Ralley  (1981). 

Since  the  critical  region  for  the  entropy  test 
involves  the  lower  (left)  tail  of  the  null  distribution 
of  H  ,  the  (asymptotic)  cutoff  point,  say  for  a  =  .05, 

IB 

corresponds  to  the  5%  point  of  the  asymptotic  null 
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distribution  and  is  seen  to  be  to  the  left  of  most  points 
in  Figure  1.  Thus,  it  is  clear  from  the  plots  in  Figure  1 
that  all  generators  under  study  except  RNCG  pass  this 
test  of  uniformity  for  all  the  values  of  m  considered. 

(Of  course  RN5  fails  for  m  <  5.) 
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3.  CONVERGENCE  OF  H  „  TO  ASYMPTOTIC  NULL  DISTRIBUTIONS 
_ m,n _ 

Dudewicz  and  van  der  Meulen  (1979)  reported  Monte 
Carlo  percentage  points  of  the  null  distribution  of  n 
for  10,000  replications  of  sizes  n  =  10(10)80(50)100  each 
using  the  generator  UNI  (using  ISEED  =  S24287  and  JSEED  = 
654345465).  A  pilot  study  showed  that  extending  this  to 
10,000  replications  of  size  10,000  each  would  be  feasible 
but  expensive. 

In  the  same  report ,  the  above  authors  also  studied 
the  normal  approximation  to  the  percentage  point  Ha*(m,n) 
given  by  (1.11)  for  a  =  .05,  and  a  variety  of  m,n  up  to 
n  _<  100.  It  was  found  that  the  asymptotic  cutoffs  deviate 
by  a  large  amount  from  the  simulated  cutoff  points,  at 
least  for  n  <_  100. 

In  this  report  we  continue  the  investigations  of  the 
above  authors  and  investigate  the  normal  approximation  to 
percentage  points  for  n  =  10,000.  For  this  we  first 
obtain  estimates  of  the  true  percentage  points  based  on  a 
Monte  Carlo  study  using  1,000  samples  of  size  10,000.  The 
resulting  Monte  Carlo  estimates  are  given  in  Table  1,  for 
various  levels  of  a(a  -  .40,  .30,  .20,  .10,  .05,  .01)  and 
m  =  1(1)10(5)20(10)40.  We  also  obtained  the  corresponding 
asymptotic  values  of  Ha*(m,n),  for  n  =  10,000  and  the  same 
values  of  a  and  m.  These  asymptotic  values  are  listed  in 
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Table  1,  together  with  the  Monte  Carlo  estimates,  and  it 
is  seen  that  -  for  such  a  large  sample  size  as  n  =  10,000 
the  asymptotic  cutoffs  are  -  for1  m  up  to  10  -  in  close 
agreement  with  the  simulated  ones  using  UNI,  and  based  on 
1000  replications. 

In  addition  to  the  above  we  plotted  the  asymptotic 

distributions  of  the  H  -  statistics  for  n  =  10,000  and 

m  ,n  ’ 

m  =  1(1)10(5)20(10)50  ,  250  in  Figures  2a-2p.  This 

1  /  3  6 

asymptotic  distribution  is,  for  n  -*■  m  =  o(nA  “  ),  and 
6  >  0,  given  in  an  explicit  form  by  ( cf .  1.10) 

l  ( -log (2m)  -  Y  +  R( 1 ,  2m~l),  gi^) ,  (3.1) 

where  j(y,oz)  denotes  a  normal  distribution  with  mean  p 
and  variance  a2.  In  each  figure  the  asymptotic  distribu¬ 
tion  function  of  n  is  overlayed  with  the  corresponding 
simulated  percentage  points  of  Hm  n  obtained  from  the 
Monte  Carlo  study  described  above  for  a  =  .005,  .010, 

.025,  .050,  .100,  .200,  .300,  .400.  It  is  seen  from  these 
figures  that  the  normal  approximation  to  the  simulated 
values  is  good  for  all  m-values  up  to  4,  bad  beyond  10, 
and  questionable  in  the  range  inbetween.  The  goodness 
depends  on  the  use  to  which  the  approximation  is  to  be 
put.  For  the  use  in  Section  5  below  it  is  bad  in  the 
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range  m  >  4.  (From  the  restriction  on  m  one  expects  good¬ 
ness  at  best  for  m  <  ^  22.) 


FIGURES 

2a-2p 

The  fact  that  the  normal  approximation  to  the  distri¬ 
bution  of  H  deteriorates  for  n  =  10,000  as  m  >  10,  as 
m  ,n 

observed  from  Figures  2a-2p,  also  has  a  bearing  on  the 
study  reported  in  Section  2.  There,  for  one  sample  of 
size  n  =  10,000  the  percentile  ranks  of  the  observed  values 
of  Hm  n  with  respect  to  the  asymptotic  distributions  are 
plotted  for  nine  random  number  generators  as  functions  of 
m.  Strictly  speaking,  this  approximation  is  not  justified 
for  values  of  m  >  ^  2  2  on  the  basis  of  Theorem  1.8. 

The  figures  2a-2p  indicate  that  the  asymptotic  distribution 
moves  more  and  more  to  the  right  of  the  percentile  points 
obtained  through  simulation  as  m  increases,  m  >_  15,  i.e. 
the  asymptotic  distribution  is  more  and  more  significantly 
stochastically  larger  than  the  actual  (simulated)  distri¬ 
bution.  From  this  it  can  be  surmised  that  for  m  >_  15  the 
actual  percentile  ranks  of  the  observed  H  -  values  in  the 
study  leading  to  Figure  1  are  larger  than  the  ones  plotted 
there,  which  are  based  on  the  asymptotic  distribution,  and 
that  the  discrepancies  get  larger  as  m  increases.  This  may 
explain  the  general  upward  trend  in  Figure  1  of  the  graphs 


12 


as  m  increases  beyond  15;  the  actual  percentile  ranks  are 
larger  than  the  plotted  ones  for  m  >  15,  which,  when 
corrected  for,  would  yield  a  more  horizontal  behavior  of 
the  percentile  rank  curves  for  m  >_  15  for  the  random 
number  generators  under  consideration. 


4. 


CONVERGENCE  OF  H 


TO  ASYMPTOTIC  ALTERNATIVE 


m  ,n 


DISTRIBUTIONS 

Expression  (1.13)  for  the  asymptotic  power  of  the 
test  against  alternatives  which  are  bounded  positive  step 
functions  on  [0,1]  was  used  by  Dudewicz  and  van  der 
Meulen  (1979)  to  compute  an  approximation  to  the  power  for 
nine  alternative  dist irbutions  for  n  =  20,  m=  1,  2,  4, 

9  and  a  =  .05.  Some  of  these  alternative  distributions 
belong  to  the  class  ?  of  bounded  positive  step  functions, 
others  do  not.  As  a  general  conclusion,  the  above  authors 
found  that  the  asymptotic  power  approximation  agreed 
poorly  with  the  Monte  Carlo  power  evaluations  for  the 
values  of  n,  m,  and  a  considered.  These  investigations 
were  continued  for  n  =  100  and  it  was  found  that  the 
normal  approximation  is  still  not  accurate. 

In  this  paper  we  extend  the  evaluation  of  the  normal 
approximation  to  power  for  a  particular  alternative  (called 
F  in  Dudewicz  and  van  der  Meulen  (1979))  for  n  =  100, 

1000,  and  10,000  and  a  =  .05.  This  alternative  which, 
belongs  to  ?  ,  is  defined  by  the  following  distribution 
function : 
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F:  F(x)  =  x 

F(x)  =  6(1/4- e/2)+x(l-6) , 
F(x)  =  x-e5, 

F(x)  =  -6( 3/4+e/2)+x(l+6) , 


if  0  <  x  <  l/4-e/2 , 
if  l/4-e/2  <  x  <  1/4  + 
if  1/4  +  e/2  <  x  <  3/4 
if  3/4-e/2  £  x  £  3/4  + 
if  3/4  +  e/2  <  x  <  1 , 


F(x)  =  x,  __ 

(with  0  <  e  <  0.5,  0  <  6  £  1;  we  used  e  =  .02,  6  =  0.1). 
In  order  to  carry  out  the  asymptotic  power  evaluation 


for  alternative  F  we  first  need  the  cutoff  points  HQ*(m,n) 


of  our  test  procedure  (1.7)  for  various  values  of  a,  m,  n, 
and,  moreover,  Monte  Carlo  estimates  of  the  power  to  compare 
the  asymptotic  approximation  with.  Monte  Carlo  estimates 
of  Ha*(m,n)  for  n  =  100  and  various  values  of  a  and  m  are 
given  in  Dudewicz  and  van  der  Meulen  (1979).  Monte  Carlo 
estimates  of  Ha*(m,n)  for  n  =  1000  (based  on  N  =  10,000 
replications,  and  for  m  =  1(1)10(5)20(10)50,  250)  and 
for  n  =  10,000  (based  on  N  "  1000  replications,  and  for 
m  =  1(1)10(5)20(10)50,  250,  500,  2500)  are  given 
in  Table  2  and  Table  3  respectively  for  a  =  .005,  .01, 


.025,  .05,  .10,  .20,  .30,  .40. 


Table  2 


Table  3 


Monte  Carlo  estimates  for  the  power  of  the  entropy- 
based  test  against  alternative  F  for  n  =  100,  1000,  and 


e/2  , 

-  e/2, 
e/2, 
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10,000  (based  on  10,000,  10,000,  and  1000  replications 
respectively)  are  given  in  Table  4  for  a  =  .05  and 
m  =  1(1)10(5)20(.l0)40  ,  together  with  the  asymptotic  power 
calculated  in  each  case  using  formula  (1.13).  Both  the 
Monte  Carlo  estimates  ("simulated  power")  and  the 
asymptotic  power  are  obtained  using  the  Monte  Carlo 
estimates  of  Ha*(m,n)  ("simulated  cutoffs")  described 
above . 

In  calculating  the  asymptotic  power  note  that  -  for 

the  density  ^  of  F  -  H(f1)  =  -.000002  ,  Var  f  log-  f^X) 

1  “ 

-.000003842,  and  y  =  .5772... 

It  is  seen  from  Table  4  that  there  is  reasonable 
agreement  for  n  =  10,000  (and  1  <  m  <  3),  but  for  n  =  100 
and  n  =  1,000  there  seems  to  be  poor  agreement  overall, 
except  for  m  =  1  at  n  =  1,000. 

Table  4 
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5. 


EVALUATION  OF  NINE  RANDOM  NUMBER  GENERATORS  BASED 


ON  A  x2~TEST 

In  Section  2  we  presented  a  one-sample  evaluation  of 

nine  random  number  generators  based  on  the  percentile 

ranks  of  the  observed  H  -values  within  the  asymptotic 

m  ,n 

distribution. 

2 

As  a  more  sensitive  approach  we  ran  a  x  -test  using 

N  =  100  sample-values  of  H  (n  =  10,000;  m  =  1(1)4) 

m,n 

from  the  same  nine  (without  RN5)  random  number  generators 

to  determine  whether  the  distributions  of  H  fit  the 

m  ,n 

corresponding  asymptotic  distributions  which,  for  m  <  4, 
have  been  verified  in  Section  3  as  being  very  close  to  the 
true  distributions.  The  results  are  given  in  Table  5. 

Table  5 
2 

For  carrying  out  the  x  -test  for  goodness  of  fit  we 
used  10  equi-probable  cells.  The  null-hypothesis  enter¬ 
tained  in  each  case  is  that  the  100  sample  values  or 
Hm  n  come  from  the  asymptotic  normal  distribution  given 
by  (3.1).  The  entries  in  Table  5  are  the  observed 
X  -values  for  the  eight  random  number  generators  and 
four  values  of  m.  Those  x  -values  which  are  significant 
at  the  .05  level  are  starred. 


o 

From  the  observed  x  -values  we  conclude  that  seven 
random  number  generators  produce  numbers  which  are  con¬ 
sistent  with  the  hypothesis  of  randomness  as  measured  by 
the  Chi-square  entropy-test ,  whereas  one  random  number 
generator  (RN3)  is  to  be  rejected  due  to  the  lack  of 

proper  distribution  of  the  H  -  statistic  on  the  basis  of 

m  ,n 

samples  of  size  n  =  10,000. 
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6.  CONCLUSIONS 

In  this  paper  we  looked  at  the  performance  of  nine 

random  number  generators  in  terms  of  the  entropy-based 

test  statistic  H  .  A  one-sample  evaluation  of  size 

m  ,n 

n  =  10,000  indicated  that  all  (except  two)  random  number 
generators  easily  pass  the  entropy -based  test,  but  it 
also  showed  that  among  the  various  generators  some  are 
closer  to  the  uniform  distribution  than  others  as 
measured  by  the  value  of  n.  One  random  number  generator 
was  excluded  from  further  study  since  its  small  period 
leads  to  rejection  at  every  level. 


The  study  of  the  convergence  of  Hm  n  to  the  asymptotic 

null  distribution  (initiated  in  Dudewicz  and  van  der 

Meulen  (197S)  for  sample  sizes  up  to  n  =  100,  and  continued 

here  for  n  =  10,000)  revealed  that  for  n  =  10,000  and 

m  10  the  normal  approximation  to  the  percentage  points 

(and  thus  to  the  distribution  of  H  )  is  good.  This 

m  ,n  ° 

comparison  was  carried  out  in  two  ways:  i)  by  a 
comparison  of  simulated  cutoff  points  with  those  found 
through  the  asymptotic  formula  (1.11),  and  (ii)  by  a  plot 
of  the  asymptotic  normal  distribution  of  H  with  the 


simulated  values  drawn  into  it.  The  latter  study  confirmed 
the  first  one. 

Next  we  investigated  the  convergence  of  H  ^  to  the 

m,n 

asymptotic  alternative  distribution  under  a  specific 

alternative  F  which,  because  it  is  so  close  to  uniform, 

is  hard  to  detect.  We  observed  that  the  simulated  power 

of  the  entropy-based  test  is  small  for  n  =  10,000  but  for 

some  m's  definitely  larger  than  the  significance  level. 

We  noticed  that  the  asymptotic  approximation  to  power  for 

n  =  10,000  as  given  by  formula  (1.13)  has  validity  only  for 

m  5  3.  In  the  course  of  these  investigations  we  obtained 

Monte  Carlo  estimates  of  the  percentage  points  for  n  =  1,000 

and  n  =  10,000  which  supplement  the  tables  already  provided 

in  Dudewicz  and  van  der  Meulen  (1979). 

Finally  we  evaluated  the  eight  remaining  random  number 

generators  on  basis  of  a  x  -test.  The  purpose  of  this  test 

was  to  investigate  whether  the  true  distribution  of  H 

°  m  ,n 

fitted  the  asymptotic  null  distributions  closely,  which  we 
know  from  the  results  in  Section  3  they  should  for  m  =  1, 

2,  3,  4  if  the  generators  are  truely  random.  On  basis  of 
this  study  the  random  number  generator  RN3  was  rejected. 
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Comparison  of  9  random  number  generators  using  percentiie  ranks  of  the 
asymptotic  distribution  of  H  for  n=10,000  and  m=l ( 1) 10 { 5) 20 (10) 40. 


Figures  2a- 2p.  The  asymptotic  distribution  function  of  for  n=10,000  and 

ffl=l(l) 10 <5) 20(10) 50, 250;  *  represents  a  Monte  Carlo  value  at 
a*=. 005,. 01,  .025,  -05,-10,  .20,  .30,  .40. 

Figure  2a. 

N=  10000 

M=  1 
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Figure  2f. 
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Figure  2p. 
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ASYMPTOTIC  POWER  APPROXIMATION  EVALUATION 
(ALTERNATIVE  :  F,  ALPHA  =  .05) 
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Table  5 


m 

UNI 

KERAND 

RN2 

RN4 

RNCG 

RN1 

SRAND 

RN3 

1 

6.0000 

11.00 

16.40 

12 .20 

7.40 

2 

4. 8000 

6 .40 

10.20 

10.60 

10.60 

3 

5. 6000 

6.00 

6.60 

8 . 00 

11.60 

4 

14.4000 

8.00 

7.60 

10.40 

19.20* 
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